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ABSTRACT 

A  model  for  an  ideally  polarizable  interface  is  proposed  where  the 
two  sides  of  the  interface  are  described  by  two  interacting  classical 
one  component  plasmas  of  different  neutralizing  background  densities 

separated  by  an  impermeable  membrane..The  statistical  mechanics  of  this 

n 

system  is  solved  exactly  in  two  dimensions  at  the  reduced  temperature  2. 

The  one  and  two  body  distribution  functions, the  potential  drop, 
the  differential  capacity  and  the  interfacial  tension  are  computed. 

The  main  feature  of  this  model  is  that  the  potential  drop  and  not  the 
surface  charge  appears  as  the  natural  external  variable. Several  sum  rules 
are  discussed, such  as  the  screening  theorems  for  multipoles, the  contact 
theorem,  and  Lippmann's  elation. 

A  general  thermodynamic  argument  shows  that  these  sum  rules  remain 

valid  in  three  dimensions  and  for  any  value  of  the  coupling  parameter. 

/\ 


I  INTRODUCTION 

In  recent  years  there  has  been  much  interest  in  the  statistical  mechanics 
of  inhomogeneous  charged  systems  such  as  the  metal  electrolyte  interface. 

Since  the  real  interface  is  a  very  complex  system, most  of  the  effort 
has  been  directed  to  understanding  models  in  which  the  metal  side  is 
an  ideally  smooth, hard, charged  surface, with  or  without  metal  boundary 
conditions  (image  forces) .Accordingly  ,  the  interface  is  not  treated  as 
an  wholy  integrated  system, in  which  the  metal  and  the  electrolyte  are 
both  treated  at  the  microscopic  level, and  this  is  surely  a  serious 
shortcoming  of  these  models.  [1] 

For  instance, if  one  considers  the  interface  between  a  liquid  metal 
(3g  or  Ga)  and  an  aqueous  electrolyte  solution  (e.g.  XC1  )  , it  appears 
that  the  main  contribution  to  the  interfacial  tension  is  coming  from 
the  interaction  between  the  metal  particles. Since, from  Lippmann's 
equation,  the  differential  capacity  of  the  interface  is  the  second  derivative 
of  this  surface  tension  with  respect  to  the  applied  potential, it  is 
very  likely  that  a  sizeable  part  of  this  capacity  is  also  due  to  the 
response  of  the  metal. 

Indeed  this  effect  has  been  estimated  recently  by  a  model  in  which 
the  metal  is  a  quantum  mechanical  jel'l,ium  [2, 3], and  more  recently 
a  jellium  interface  with  an  ionrdipole  mixture  has  been  consider ed[ 4, 5]  . 
However,  this  last  model, albeit  realistic,  is  very  difficult  to  handle  and 
we  are  far  frcm  an  accurate  numerical  solution  of  the  statistical  mechanics 
of  this  modal. we  must  remark  that  even  computer  simulations  cannot  be 
performed  for  such  a  complex  case  because  of  the  quantum  nature  of 
the  metal, and  because  the  discrete  solvent  of  the  ionic  solution 
will  require  a  number  of  particles  that  is  too  large  even  for  modern 
computers. 


z 


k  It  seems  therefore  interesting  to  study  simpler  models  that  may 
not  be  as  realistic  #  but  still  treat  the  whole  interface 

at  the  microscopic  level  #and  are  more  tractable  mathematically .In  particular 
we  would  like  to  mimic  two  different  situations  that  occur  in  electrochemistry 

a)  the  non  polarizable  interface# in  which  the  potential  drop  across 
the  interface  is  given  by  Nernst's  equation. The  classical  example 

of  this  case  is  the  calomel  electrode  and  its  behaviour  is  explained  by 
the  continuous  exchange  of  charged  particles  between  the  two  sides  of 
the  interface. 

b)  the  polarizable  interface# in  which  the  potential 

drop  does  depend  on  the  charge  distribution  at  the  interface. 

The  classical  example  of  thi?  case  is  the  mercury-FNa(ag)  interface. 

# 

For  low  values  of  the  potential  drop  across  the  interface 
a  change  in  the  externally  applied  potential  will  not  simply  produce 
an  increase  in  the  ohmic  current  intensity# but  it  will  polarize  the 
surface# producing  an  excess  surface  charge  on  the  metal  and  in  the  solution. 
Above  a  certain  value  of  the  potential  however  chemical  reactions  occur  and 
electrons  flow  across  the  interface:  At  this  point  the  polarizable  electrode 
becomes  urroolarizable  ,and  behaves  much  like  case  (a). 

In  the  oresent  work  we  consider  a  2-dimensional  model  in  which  the 
two  sides  of  the  interface  are  represented  by  classical#one  component  plasmas 
(OCP)  with  different  background  densities. Of  course  this  is  not  a  realistic 
representation  of  the  metal  electrolyte  inter f ace #but  this  model  has  the 
feature  of  being  exactly# analytically  solvable  :Jancbvici  and  coworkers [6,7,8] 
have  shown  that  at  the  special  value  of  the  plasma  coupling  parameter 
(or  reduced  temperature  )  -2, the  partition  function  and  all 


th«  correlation  functions  of  the  2- dimensional  OCP  can  be  obtained  exactly 

On  the  other  hand, as  ve  shall  see  below, this  model  reproduces  some 
of  the  main  features  of  real  interfaces  and  may  be  useful  in  providing 
insight  about  exact  relations  such  as  sum  rules  and  the  subtle  screening 
behaviour  of  Coulomb  systems  [9]. 

Recently  a  solution  was  found  for  such  a  model  of  the  interface 
where  the  charged  particles  (ions)  are  free  to  move  everywhere 
(10 ,11]. In  this  case  the  potential  drop  across  the 
interface  depends  only  on  the  difference  in  the  chemical 
potentials  of  the  bulk  phases,  and  any  surface  charge  distribution  is 
perfectly  screened  by  the  mobile  ions, giving  no  contribution  to  the 
potential  drop. Clearly,  this  model  mimics  the  behaviour  of  the  non 
polarizable  electrode. 

In  the  present  paper  we  present  a  generalization  of  this  model 

0 

which  mimics  the  behaviour  of  the  polarizable  interface.  Por  that  purpose 
we  introduce  an  impermeable  membrane  between  the  two  media, which  prevents  the 
ions  from  moving  across  the  inter face, while  allowing  the  electrostatic 
interactions  among  all  ions  and  background  charge. Indeed  a  completely  • 
impermeable  membrane  is  unlikely  to  exist  in  nature,  and  there  is  always  some 

small  but  finite  current  which  leaks  across  the  interface. In  a  limited  range 

/ 

of  applied  potential, however,  a  system  like  Hg-PNa(aq.) ,caa  be  considered 
as  ideally  polarizable, because  the  leakage  current  is  extremely  small. 

We  must  point  out  that  in  real  systems,  the  concept  of  an  .ideally 
polarizable  electrode  does  not  imply  the  existence  of  a  zero  width  , 
mathematical  surface  ,but  there  could  be  a  region  of  limited  extent  where 
all  particles  could  coexist  (12] . 

We  shall  not  consider  this  case  here  because  it  . is  more  complex  than 
the  former  one. 


In  section  2  we  describe  the  model  and  the  method  of  solution 
Section  3  is  devoted  to  the  discussion  of  1  and  2  particle  density 


functions, and  related  sum  rules. In  Section  4  we  compute  the 
total  potential  drop  across  the  interface,  and  show  that  in  our  model 
this  is  indeed  the  natural  external  variable, while  the  surface 
charge  is  calculated  by  means  of  a  closed  formula  from  it. In  this 
sense  our  model  differs  from  most  of  the  literature,  where  the  charge 

i 

is  the  externally  fixed  variable  and  the  potential  drop  is  calculated 
from  it. The  thermodynamics  of  the  system  is  discussed  in  Section  5  and 
pertinent  sum  rules  are  given: in  particular  Lippmann's  equation 
and  the  contact  theorem  for  this  case  are  deduced.  The  last  section 
is  devoted  to  a  brief  discussion  of  the  results. 

A  short  communication  of  these  results  has  been  presented  elsewhere [13] . 

II-  MODEL  AND  METHOD  OP  SOLUTION 

Following  the  method  of  Jancovici  [8]  we  start  with  a  system 
of  circular  symmetry. In  our  case  it  consists  of  two  concentric  circles 
of  radii  R^  and  Rt 

The  inner  region  (0<r<^,  )  has  background  charge  density  -ec , where 
e  is  the  elementary  charge, and  contains  N.  particules  of  charge  e  . 

This  circle  is  also  the  impermeable  surface  that  prevents  the 
charged  particles  from  changing  sides. 

We  want  to  allow  for  an  excess  charge  density  -ec'^  along  the  wall 
then  we  must  have 


N.  -  ^ 


(2.1) 


The  outer  region  (*H  <*<%.)  has  background  density  -eot^  f excess 
charge  density  eo5£.  and  contains  particules  of  charge  e.'lence 


hft  x  +  *  O'**  (2.2) 

the  system  is  thus  globally  neutral. The  thermodynamic  limit  is  obtained 
letting  the  number  of  particles  and  the  radii  go  to  infinity  with  ^ 


(2.2) 


and  0"  constant . 


The  plasma  coupling  parameter  is  f*  &  A  where  *l/kT  is  the  thermal 


Boltzmann  factor  and  T  is  the  temperature 


The  Coulomb  potential  between  two  particles  at  a  distance  r  is 
O'Cr'l  -  _  where  L  is  seme  scaling  length. 

The  first  step  is  to  compute  the  total  potential  energy  of 
the  system.  We  find  ^  U 

v  -  v.  *  **  l-  r  *  s  t  + 

*  CV-.V*?  Z  (*.*) 

(•ft/,*!  *  J 

where  'n  is  the  position  of  ion  i,N«H,+^and  V#is  a  constant  background  term 

*  «i  ( r\  -  \  ca.k\ 

Note  that  the  excess  charge  density  ecr  does  not  appear  explicitly 
in  (2.3)  and  (2.4). 

In  computing  the  canonical  partition  function,  the  angular  integrations 
can  be  readily  performed  expanding  a  Vandermonde  determinant  17,8] 
and  we  get 


?  =  —  ,  {&•+* 

*  WJ  K 

.1  •  I  " 


• t <? 

CP 


where 


<?  is  the  permutation  operator  of  all  particles. 
This  can  be  written 


tf-t 


*  «*vrz:  nu  /i 

£P'  n.i  •  %  V*^ 


where  r»ts  n.  + 

and  is  the  permutation  operator  that  exchanges  the  J  hels  {0,^-1} 

with  labels  {1^  ,N-1}  in  all  possible  ways. 

Introducing  the  functions 

M.  9* 


and 


(5-T'l 


[  y  (V  .*»*'«  1  -  y(v'*SR^] 


where 


rMA*)*  e'V J* 

is  the  incomplete  gamma  function, we  get 

=  e*^#,n:rf  CT  [  v£>V>  v£  Vi . .  4>V»yj 


H.-i 


CP’  '  '  as 

This  can  be  written  in  a  more  convenient  form  factorizing  the  term 


IT  vi>V^  TT  <4>Vv  . 


tfe  find 


where 


Z^  r  Z*  rE 

M  1 - 


(J.O 


and 


2^  s  e  tc^  JT  /  /  c±>V»c\ 

n*«/.  ' 


(*.  U>\ 


is  the  generating  function 
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NU 


I - I  ”  -1 


£ _  /  jc 

v4  m.s  V  xTm 

M,-i  M-2 


+  i _  / _  jCKvXM1 

n.*'*/V"-»  XC^Ay^rv' 


where 


/ _  l_ _ 

ni,H/Vn.,Vfl»  »v  l,itl;V’'/V",t  txcM  xov^Jtev^ 


X<r>^*  +W4>V* 

In  (2.11)*  it  appears  that  in  the  limit  R^  ,R^  ->©©,  the  dominant  values 

/ 

of  the  arguments  of  the  summations  are  those  with  indices  close  to  N  . 

The  incomplete  gamma  function  can  be  expanded  in  this  limit  using  the 
asymptotic  formula  [6] 

r  ('•♦•/ *1*? ^  «  »■  [  $-ii) 


(5. 'A 


where 


t  = 


and  s  i  J  <f”  Ju. 

is  the  complementary  error  function. 

This  yields, after  some  manipulations 


aw=  xcy\  -r 

j(  (t*\  =  ra  <2  &rfe  (m^ 


where 


(2^ 


and  ms  VJ{u_  ° 

I  g,  ^ 

The  deviations  of  ul-i  from  unity  are  due  to  surface  contributions. 
In  a  first  approximation, we  can  simply  neglect  these  contributions. 

This  will  be  a  mean  field  approximation  (UFA)  .The  neglect  of  the 

rn 

variation  of  the  indices  in  the  computation  of  i_j  constitutes  a  second 
approximation  which  will  be  found  to  be  surprisingly  accurate. 


o 


a)  The  mean  field  approximation  (HFA) 
Since  S  »1  we  have 


W— t 


that 


?US  2#  s  TC*  ll  I  l 

**  nso  ft*  */, 

This  approximation  has  a  physical  meaning: It  is  clear  that 

2°  is  the  partition  function  of  the  uncoupled  system  where  one 

neglects  the  correlations  of  the  charge  fluctuations  between  the  two 

sides  of  the  inter  face. In  this  sense  it  is  a  bona  fide  mean  field 

approximation. The  system  is  the  superposition  of  two  subsystems :The  inner 

OCP  enclosed  by  a  ring  of  radius  R^with  charge  density  £0/^ 

and  the  outer  OCP  enclosing  the  charged  ring  of  radius  and  charge 

density  —  . 

These  cases  have  been  already  solved  in  the  literature!  3 ,1-1]  ,so 
that  we  merely  quote  the  results  for  the  partition  function  and  the 
one  body  density. For  the  inner  circle  subsystem 

\  -  ih  **  =  ^  ft  i>\A  «.*  \ 

^  nap* 


where 


is  the  constant  background  term 


V‘°  :  7  [  ^  ] 

and  for  the  outer  ring  subsystem 


2  -  -f 

*"*r  —  — — ■ 


where 


±- 

V  J 


Tf 

»» fC 


(MO 


(*•  I?'! 


A 

t  «,**,'**.  * 

and  the  particles  have  been  relabelled  from  {0,Ht-l}  to  {K^W-l}  . 

It  can  be  easily  verified  that  ...  a/./o  ,©\ 

-fit*  -pvYi  *v,  t 

?'  .  _e_j - i  CMO 


*T% 


I 


hence 


In  the  limit  of  a^olanar*  interface  (R.,R.->oo),  the  corresponding  one 

£  ■* 


particle  density  profile  is  [8,14] 


1U,  J  iYl  jy. 

**  r4k. 

.  .  f  ftftcfiC-fr  j  , 

Id 

with  origin  of  coordinates  at  the  wall .  (x-r-R^ )  . 


X  *  o 


X  *  o 


C*2l'\ 


In  two  dimensions 


>2  corresponds  a  somewhat  weak  coupling 


for  the  bulk  ions, so  we  can  expext  that  this  mean  field  approximation 

will  yield  sensible  results  for  not  too  large  surface  charges. 

We  think  that  for  any  realistic  model  of  the  interface  this  will  be 

«  // 

a  first  approximation. If  we  also  neglect  the  contributions  of  the  metal  side 
(say,  the  RHS  of  the  interface),  then  we  will  recover  a  model  of  the  interface 

^  it 

similar  to  the  one  discussed  in  most  of  the  literature (the  primitive  electrode 
b)  Secojnd  approximation. 

A  better , perhaps  not  as  transparent  physically , approximation 
can  be  obtained  by  neglecting  end  terms  in  the  summations  (2.11) . 

Using  the  identity  ^ 

94  94  -  J 


ZL  .'LL 


can  be  written 


*9  L  L  ™ 

*4  x-t  ,  Pf-I  ,  tf 


Xt  ,  n,-H  V-4  ,  »/ 

L  J  XM 


*2  >  +  ....  (4.ti) 

•>,**,  nt* */***•'  *»,♦-»  •5CCwj') 

which  shows  that  the  terms  neglected  by  extending  the  sums  at  the  ends 


of  the  intervals  are  small  (see(2. 13) ) .Then  we  write 


w, 


X<r.\ 


Z  ,  4 .  an  —  *  2_  z 

nt& i  OCC*r\  %*1  / *  v^iW  x^^i\ 

►k  ,  ft,.!,  »t.t  _ _ 

^  ^  J  *_..  (j.t^ 

In  the  limit  N,  ->oo  ,  the  sums  can  be  replaced  by  integrals  over 

the  new  variable  -  i«  and  using  the  formula 

\fHt 

.AA-A  - 


we  get 

r — i 

i _ i 

*£  +  lK| 

+  •  •  •  e 

<WX 

C*0r 

where 

r»fA i,  , 

/  -iiL 

GO 

And 

«,"  / 

(■?.  u\ 


which  can  be  recognized  as  the  expansion  of  the  modified  Bessel  function 

xo  U'R  ^>0  [15]  . 

Asymptotically 
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We  shall  discuss  this  interesting  approximation  below, 
c)  Exact^  solution.  . 

Instead  of  using  (2.11)  to  compute  2  the  exact  partition  function 
of  the  coupled  system, we  could  represent  it  in  a  much  more  compact  form 
using  a  device  similar  to  a  grand  canonical  ensemble  calculation.Consider 


now  that  all  particles  are  indistinguishable  and  can  move  freely 
throughout  the  entire  system, but  assume  that  the  particles  of  the  inner 
region  have  a  fugacity  different  from  that  of  the  outer  region, which 

we  take  as  the  reference  and  set  it  to  the  value  of  1. 

It  is  straightforward  to  see  that  the  partition  function  2^  is  the 
coefficient  of  4  in  the  polynomial 

W-/ 

^  °  TT  i  yf'c*')  ♦  (?.**\ 

which  can  also  be  written 


Using  the  theorem  of  residues,  we  finally  get 
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which  can  be  evaluated  asymptotically  by  the  method  of  steepest  descents. 
The  path  of  integration  is  a  circle  centered  at  the  origin. 
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The  saddle  point  is  located  at  the  positive  real  ^axis  .The  location 
of  the  saddle  point  is  found  by  solving  the  equation 
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Taking  again  the  Unit  of  an  infinitely  large  system, and  changing  to 

the  variable  t  «  *i ^t*-" *- 

V  2.*^, 

this  equation  becomes 
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For  m»l  we  have  ^)C(rfc^  s  OOC^""*  ' 

so  that 

arN3  the  general  case  this  equation  must  be  solved  numerically 

for  given  values  of  the  surface  charge  density  ee*. 

Now  we  take  -  ^ejrp/t0^  and  then  integral  (2.32)  yields 
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Now  we  take 
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Equations  (2.9)  ,(2.35)  and (2. 36)  are  the  exact  answers  for  the 


partition  function  in  the  fully  coupled  case. 

fts  one  would  expect  in  the  asymptotic  limit  of  a  large  system,  -fit EL 
is  proportional  to  the*surface*  of  the  inter face. Hence  the  Helmholtz 
free  energy  of  the  coupled  system  differs  only  by  a  surface  term  kf* frctn  the 
free  energy  of  the  uncoupled  system: 

&f*  s-  Vr-fU. 
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This  equation  gives  the  thermodynamic  properties  of  our  model:It 

% 

will  be  discussed  in  detail  in  section  5. 

A  very  good  approximation  of  the  exact  solution  of  the  saddle  point 
equation  (2.34)  can  be  obtained  if  we  assume  that 
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and  expand  in  the  small  variable 


(2.34)  .We  find  to  first  order 


3.  *  0»* 

where  and  Kj  are  defined  in  (2.25) .The  same  procedure  applied  to  (2.37) 
yields 


A/s  =  -  si  tfS, 


But  this  is  also  the  result  obtained  directly  from  (2.27)  .Hence  this 
approximation  is  strictly  equivalent  to  the  one  of  subsection  (2-b). 


3. DENSITY  PROFILE  AND  PAIR  CORRELATION  FUNCTIONS 


The  same  methods  used  in  obtaining  the  partition  function  can  be 
used  in  obtaining  the  one  and  two  body  correlation  functions, 
a)  DensjLt^  orofjle^ 


From  equation  (2.28)  we  find  that  the  density  profile’ is  for  r<R 
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where  A  is  the  coefficient  of 
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Hence, after  some  manipulations, we  get 
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Similarly, for  r>R 
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In  the  plane  interface  limit  we  use  again  the  asymptotic  expansion 
of  the  incomplete  gamma  function  and  replace  the  sums  by  integrals  over 
t  .The  asymptotic  behaviour  of  the  contour  integrals  is  again  obtained 
by  steepest  descents. The  saddle  point  is  located  at  the  same  position 
as  before. 

After  some  manipulations  we  obtain 
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which  can  be  written  in  a  more  compact  form 
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Although  the  physics  of  the  .present  model  is  not  quite  the  same  as  that 
of  the  permeable  membrane  [10/111, it  is  interesting,  to  note  that  we 
recover  the  density  profile  of  this  case  when  ^»l.In  this  former  case 
the  density  of  mobile  ions  was  always  continuous. In  the  present  model, however, 
there  is  a  discontinuity  in  the  density  profile  at  the  interface: 
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Pot  m-1  (o^«  o^»  c<  >  these  contact  values  can  be  computed  explicitly 
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For  m«l  and  0* «o  one  obtains  the  expected  result 
X  )  ■  constant 

which  means'  that  the  system  behaves  as  if  there  would  be  no  wall. 

* 

The  approximation  of  subsection  (2~b)  is  of  interest  be'cause  the 
density  profile  ^p(x)is  the  sum  of  the  density  of  the  MPA  plus  a  perturbation, 
Osing  (2.38)  we  find 
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wh«*«  i*  given  by  (2.’ 21). 

The  integrals  in  (3.5)  and  (3.9)  have  been  computed  numerically 

for  the  case  *^-1,  ^-4  and  O'*  — l.The  density  profiles  are  shown 

in  figure  1, together  with  that  of  the  uncoupled  system.  Indeed, approximation 
of  subsection  (2-b)  (or  3.9)  is  very  accurate. 
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The  electroneutrality  conditions  on  each  side  of  the  interface  imply: 
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In  Appendix  A  we  show  that  these  relations  can  be  derived  from 
(2. 34), the  equation  for  the  saddle  point. This  constitutes  an  alternative 
way  of  calculating  ^  from  the  suface  charge  density  cff* . 

The  approximate  solution  of  section  (2-b)  satisfies  this  condition  exactly. 


b)  Two  body  correlations. 

To  discuss  the  2-body  correlation  functions  we  need  to  consider 
the  three  cases  of  both  ions  on  side  l,both  ions  on  side  2  and  one  ion 
on  each  side  of  the  interface. S ince  the  derivations  are  very  similar, we  shall 
give  the  details  for  the  last  case  only 

It  will  be  convenient  to  use  a  complex  number  representation  of  the  position 
of  the  ions  (Janccvici  16])  2«x+iy*r  exp(  i®)  .Then  similarly  to  (3.-1  )  we  have 
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where  k  is  the  coefficient  of 
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in  the  polynomial 
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After  some  manipulations  we  get  the  following  expression  for  the 
truncated  two  body  density: 
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In  the  plane  surface  limit  , again  by  the  method  of  steepest  descents 
we  obtain 
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where  is  the  sane  saddle  point  as  before. Using  the  expression  (2.12) 

for  the  incomplete  gamma  function, and  replacing  the  sum  by  an  integral 
over  t,we  find 
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where  y*^  -y^  .In  the  general  case  we  have 
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where  g(x)  is  given  by  (3.6). 


This  result  becomes  identical  to  the  one  for  the  permeable  membrane(10, 

when  g  -l.By  an  analysis  similar  to  that  of  Jancovici  [8,10]  we  can  see 
that  in  the  exact  solution  there  will  be  no  long  ranged  correlations  along 
the  interface:  decays  exponentially  as  a  function  of  y.In  the  MFA  however 

the  pair  correlation  function  along  the  interface  will  decay  only  as 
because  in  this  approximation  each  side  of  the  interface  is  like  an  OCP 
face  an  insulating  , charged  surface  [8]. 

Our  results  are  in  accordance  with  Jancovici's  conjecture  that  short 
ranged  decay  of  the  pair  correlation  function  is  a  general  feature  of 
the  interface  between  two  conducting  media[10]. 
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The  behaviour  of^fx^  ,  x^,y  )  for  the  same  case  as  figure  1  is  shown  in 
figures  2,3  and  4  for  different  locations  of  the  particule  2  in  the  RHS  of  the 
membrane.  In  figure  4, the  correlations  are  almost  these  of  the  bulk  phase[6] 

Because  of  the  fast  decay  of  the  correlation  function  in  all  directions 
the  correlation  functions  will  obey  the  theorems  concerning  the  nonexistence 
of  multipoles  ,  that  is  any  finite  size  fixed  multipole  will  be  perfectly 
screened [16]  .In  appendix  B  we  show  explicitly  that  both  the  monopole  and 
dipole  sum  rules  are  indeed  satisfied. 


4-  POTENTIAL  DROP  AND  DIFFERENTIAL  CAPACITY 


a)  jPotential  drop 

The  potential  drop  across  the  interface 
can  be  calculated  by  computing  the  first  moment  of  the  charge  profile 
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The  integration  can  be  performed  very  easily  for  m»l,with  the  explicit 
result 


Femembering  that  •g.  was  introduced  as  the  ratio  of  the  fuga cities 
on  beth  sides  of  the  interface, it  is  evident  that  there  must  be  a  relation 
connecting  the  potential  drop  to  .Indeed  we  have 
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where  yU4  and  are  the  chemical  potentials  in  the  two  bulk  phases. For 
r*  »2  one  has  *  - 
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hence  (4.3)  becomes  a.  tm  e/**p  ^  ^  Y/*  ^  C4«•^^ 

Relations  (4.3)  or  (4.5)  reflect  the  fact  that  an  impermeable 
membrane  renders  invalid  the  internal  thermodynamic  equilibrium  condition 

*  €  C4* 

In  the  permeable  interface  case  (4.6)  holds  [17]  so  that  necessarily  £^*1. 

This  explains  the  results  found  above  for  the  one  and  two  body  density 
functions. 

It  is  now  clear  that  our  model  behaves  like  an  ideally  polarizable 

interface. While  in  the  permeable  case  the  potential  drop  always  satisfies  (4.6) 

this  is  not  true  anymore  for  the  present  model,  and  &y>can  be  considered 

as  an  external  variable:For  a  given  value  of  ^,we  compute  the  density 

profile  from  (3.4)  and  then,  the  surface  charge  density  from  the 

»  >< 

electroneutrality  condition  (3.10)  .On  the  contrary,  m  the  primitive 
electrode  (a  plasma  near  a  charged  wall)  the  surface  charge  is  the  natural 
external  var iable.This  is  an  interesting  feature  of  this  microscopic 
model  which  compares  well  with  the  experimental  situation  where  the  potential 
drop  across  the  interface  is  fixed  externally. 

The  approximation  of  section  (2-b)  gives  the  following  expression  for  ky 

*Y°  +  *  iQgf  X (r{v^  " 

where  is  the  potential  drop  of  the  MFA. 
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(4.7)  does  not  satisfy  the  exact  relation  (4.5). 

Figure  5  shows  the  relation  between  the  potential  drop  and  the 
surface  charge  for  the  primitive  electrode  and  for  the  exact  solution  of 
the  full  coupled  model. In  the  scale  of  the  figure, the  MFA , approximation  (2-b) 
and  the  exact  solution  are  indiscernable. However , the  potentials  of  zero 
charge  are  zero  for  the  MFA  ,.142  for  (2-b)  , and  .114  for  the  exact  solution. 

For  negative  charges  the  behaviour  of  the  primitive  m^'el  Is  totally 
different  from  the  coupled  model. 
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b)  Differential  capacity. 

a  quantity  of  central  interest  in  experimental  electrochemistry  is  tne 
differential  capacity .Th is  quantity  has  often  been  used  to  characterize 
the  polarizable  interface.  It  is  defined  by  C  “  ~ 

Since  from  (4.5)  1 

HE  -  -  Is  ! 
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we  can  calculate  C  as  a  function  of  ^  and  hence  as  a  function  of 
the  potential  drop  by.  Using  the  electroneutrality  relation  (3.10)  and  (3.4) 
we  find 

C  (  s  _S.£.  ifl?, 

7T 

For  m*l  and  (|  *1)  ,  the  integral  can  be  performed  explicitly 

and  gives 
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In  figure  6  we  show  the  dependence  of  C  with  the  potential. It  is  amusing 
to  note  that  the  shape  of  the  curve  has  a  hump  instead  of  the  usual  dip 
of  the  theory  of  the  diffuse  layer  [18,1?]. This  hump  is  present 
in  many  experimental  curves,  and  has  generated  considerable  controversy 
about  its  origin.The  explanation  of  the  hump  in  our  model  is  based 
on  the  fact  that  because  the  two  sides  of  the  interface  have  always 
opposite  charges,  the  mobile  ions  on  one  of  those  sides  are  pushed 
away  from  the  interface. The  distance  between  the  center  of  the  excess 
charges  of  the  positive  and  negative  sides  of  the  interface  will  increase 
with  the  charge  (or  the  potential) .The  capacity  is  therefore  maximum  when 
this  distance  is  smallest, which  occurs  near  (but  not  exactly, except  in  the 
case  m«l)  the  point  of  zero  charge  (p.z.c). 

In  the  MFA  the  capacity  curve  has  this  same  qualitative  behaviour, but 
notice  that  this  effect  does  not  exist  in  the  primitive  electrode 
(  C ->oo  for  k<f->-co) . 
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5-  THERMODYNAMICS  AMD  SUK  WJLES 


In  the  particular  geometry  of  our  model*  the  excess  Helmholtz  free 
energy  of  the  system  can  be  separated  , f or  large  ,  into  a  bulk  contribution 
which  is  proportional  to  R*  ,  and  a  surface  contribution  which  grows  as 
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In  the  MFA  the  surface  excess  free  energy  is 

•  p  ^ 

where  [14]  . 

.  /»* 

ia  4,  t  And  Dj  8 

(In  [14]  the  first  term  of  the  right  hand  side  of  (5.3)  does  not 

/ 

exist  because  Smith  considers  the  excess  free  energy  caused  Say  the 
increase  of  the  surface  charge  density  from  0  to  , while  keeping 

the  number  of  particles  constant) 

For  the  coupled  system  we  have 
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where  is  given  by  (2.37). 
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•)  j3urface  tension  and  Lippmann'a  equation 


Consider  now  the  derivative  of  the  surface  energy  with  respect  to 


the  surface  charge  density 
Prom  [14] ,  we  know  that 
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where  is  the  potential  drop  in  the  MPA,given  by  (4.8), 

For  a  charged-  fluid  in  the  presence  of  a  charged  hard  wall,  relation 
(5.5)  is  merely  a  consequence  of  the  linearity  of  the  Hamiltonian  in  the 
external  field  due  to  the  charge  density  on  the  wall  [14,20] 
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For  our  model  eq. (5.  6)  is  not  an  exact  one  because  we  are  ignoring 
fluctuation  terms  due  to  correlations  of  charge  distributions  across 
the  inter  face.  Indeed  (5.6)  is  the  MFA  Hamiltonian. 

Now, from  (2.37)  we  have 


where  we  have  used  (2.34)  to  eliminate 
Hence 
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and  using  (4.5)  we  get  finally 


(sA 


A- A  -  f 

(sa 


- 


(5-3  \ 


A  *A  + 


which  is  similar  to  (5.5)  ,but  where  the  free  energy  andYpotential  drop 
are  now  the  exact  quantities, rather  than  approximations. 

Therefore, the  fluctuation  terms  neglected  in  (5.6)  make  no  contribution 
to  this  thermodynamic  relation. We  remark  also  that  for  the  permeable 
interface, because  of  (4.6) ,eq  (5.5)  simply  means  that  the  system  adjusts 
to  minimize  its  surface  free  energy. 

As  shown  elsewhere  in  [21]  a  relation  of  the  kind  of  (5.3) 
implies  the  Lippmann  equation  of  electrochemistry  [18]  . 

Consider  the  interfacial  tension  y  of  our  system, defined  by 


0  *  'Da  J/tc/ri^a 


(s.m^ 


where  (T  is  the  volume  of  the  system, A  is  the  area  of  the  interface, and  Q 
the  surface  charge  (  Q»edA) .We  have  the  following  relation  between  y  and 

jr  =  /$  +  jfOv/M  -  it  C*-'A 

(note  that  y-  and  are  identical  only  when  0* «0)  .Hence 
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where  (5.5)  has  been  used  .We  finally  obtain, consider ing  as  our  external 


variable 
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which  is  the  standard  form  of  Lippmann' s  equation  [is] . 

We  note  that  we  have  given  here  a  statistical  mechanical, microscopic 
derivation  of  this  equation, for  this  particular  model .The  usual 

proof  of  this  equation  involves  only  thermodynamic  arguments/al though 

there  have  been  recently  several  attempts  of  microscopic  derivations  [12,22-24] 


Quite  surprisingly  it  can  be  shown  analytically , that  the 
approximate  surface  free  energy  (2.40)  also  verifies  Lippmann's  equation. 

Pigure  7  shows  the  dependence  of  the  surface  tension 
calculated  from  (5.4)  and  (5.11)  with  the  applied  potential. This 
curve  is  known  to  the  electrochemists  as  the  electrocapillarity 
curve. As  in  the  case  of  the  differential  capacity , this  curve  also 
mimics  real  systems  behaviour  with  the  characteristic  inverted  parabolic  shape 
Since  Lippmann's  equation  implies  that  the  surface  charge  is  the  slope 
of  the  electrocapillarity  curve  the  point  of  zero  charge  corresponds  to  its 
maximum. 

5-b)  Other  sum  rules 

For  an  0 CP  near  a  charged  hard  surface  the  following  sum  rule  is 
true  [21] 

-  ~€  J  [ifi*)-  fC*9)]**  (S-1*) 

where  jo  is  the  bulk  density  and  c^(x)  is  the  local  electrostatic 
potential. 

In  the  present  model  it  is  cl4ar  that  relation  (5.14)  holds  for  each 
side  of  the  interface  separately  in  the  MPA.  We  show  in  appendix  C  that 
it  is  also  true  for  the  exact  solution. Bence  we  have 
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where  we  have  used  Poisson's  equation. These  are  new  sum  rules  for  this  model 


involving  the  second  moment  of  the  charge  density  distribution 


49 


r - 

These  relations  are  also  valid  for  the  case  of  the  permeable  interface* 

When  m*l  {*,»  •<  •  <  )  we  have  the  simple  relation 

<eo  w 

A  last  sum  rule  is  the  so  called  'contact  theorem^.Por  the  permeable 
interface  it  has  been  shown  that  [17] 

s  1£«  C<fC°)~  »f  Cf*®')]  fS.  17^ 

where  P.  and  P.  are  the  bulk  pressures  on  both  sides  of  the  interface. 

The  presence  of  the  impenetrable  membrane  modifies  this  relation. We  now  have 

as  shown  in  appendix  C.Note  that  the  surface  charge  density  does  not  appear 
explicitly. When  m«l  (5.1$)  gives 

~  ~  LY/e  *  % 

which  is  equivalent  to  (3.8)  the  equation  for  the  jump  in  ion  density 
at  the  interface. 

After  some  straightforward,  but  lenghty  manipulations, it  can  be 
shown  that  the  approximate  solution  (2-b)  has  the  remarkable  feature 
of  satisfying  both  (5. IS)  and  (5.18). 

5-c)  General  derivation  of  the  preceding  sum  rules.  . 

The  sum  rules  of  the  preceding  section  are  also  true  for  the 
general  , three  dimensional  interface  of  two  OCP  at  an  arbitrary  coupling  . 

This  can  be  shown  by  an  extension  of  the  argument  used  elsewhere  [25] 


! 


which  consists  in  studying  the  variation  Sf  of  the  Belmholtz  free  energy 
produced  by  a  variation  of  the  background  density  TtlMor  of  some 
external  one  body  potential  v*p(j?j,while  holding  the  number  of  particles 
and  the  temperature  fixed. 

Proa  the  definition  of  the  Helmholtz  free  energy  we  have  in  general 

SF  s  _/ <T ^**cr  )  «.  y>cF)  JViCr^WF  (S.to\ 

The  second  term  of  the  RHS  of  (5.20)  is  due  to  the  constant  neutralizing 
background. 

Consider  now  a  system  consisting  of  two  slabs, one  of  background 
density  !*•  extending  from  -L<x<0,and  another  slab  of  density  going  from 

ft  It 

0<x<L. Eventually  we  will  take  the  limit  L->oo.The  impenetrable  membrane 
is  located  at  x»C  and  can  be  represented  by  the  external  potential 


«*ri-  * 

^  *  ' 

o 


'V.  >  o 


•x.  <  o 


for  all  particles  of  the  left  hand  side 


f  o  x  > 

cA  --  \ 

*  (  •+  »0  yt  < 


for  all  particles  of  the  right  hand  side 


The  electroneutrality  conditions  imply  that 


M,: 

1  1L 

K  -  o-) 


where  A  is  the  area  of  the  interface. 


Consider  now  the  variation  JFin  the  free  energy  when  the  position  of 
the  impenetrable  wall  is  displaced  by  an  infinitesimal  distance  $L 


while  holding  tr*  and  T  at  a  fixed  value. 

The  left  hand  side  now  goes  from  -L  to  Jl  and  has  background  density  ^(*1  * 

The  right  hand  side  extends  from  L  to  L  and  has  density  -L/V-  r<\ 

7C  4  C 


Since  A,N  ,N  ,  and  0“  are  constant, we  have 
^  z 

-t  S'*!  1.  =.  S L  —  L  s. 


(«.z*  'N 


and  the  volumes  of  each  side  of  the  interface  change  by  the  amount 


£<r  =.  -  r<j”  *  a  s*l 


Cs.u\ 


(5.20)  becomes 


o  oo 

SF  =  ft  [  j  p(*1  f  J j>Q>\  ] 

°  “*  (S.iJ  ^ 

Wcltln9  f  ‘'-pf-/5' f/’w I  *  *  ° 

\  ^  c«-“\ 

L  **>-p  [-  <ft  c«\  ]  x  >  o 


X  >  o 


and  using  (5.2.^) ,  we  find  after  some  calculations 


JF:  -  AtT  [  -  JH^]  ft.  -  «  A  Si:  |  ~ 

<1 

Introducing  the  functions  y£(x)  is  a  standard  device  in  the  theory  of  fluids 
with  hard  repulsive  cores. The  essential  feature  of  these  functions  is  that 
they  must  be  continuous  across  the  interface. 

In  equation  (5.27)  we  see  that  the  free  energy  variation  has  a 
volume  and  a  surface  contr ibution, which  are  readily  identified  with 
thermodynamic  quantities: 

£f  «  -  7;  TtT  _  T \  *  A  5/ 


(S. 


The  presence  of  the  surface  term  in  (5.28!  is  again  a  manifestation  of  the 
existence  of  a  neutralizing  background. A  similar  terra  will  also  be  present 
if  one  of  the  sides  is  a  quantum  mechanical  jellium  [26] 

From  (5.23)  and  (5.24)  we  get  for  (5.28) 

Sf  =  -  ft  ft  -  Vi  -  A  f  2£‘  - 

Comparing  the  volume  term  of  (5.27)  and  (5.29)  we  get  the  contact  theorem 
(5.18).  The  surface  term  yields 

t  s  o 

^  ~  ^  «,  /  li'  J [fQ*' \  - 

1  .  •>  o 

(,‘S.lo\ 

This  relation  is  weaker  than  (5.15)  .We  conjecture,  however,  that  relations 
(5.15)  and  (5.16)  are  valid  for  all  couplings  in  two  and  three  dimensions. 

6 -CONCLUSIONS 

In  the  present  paper  we  have  solved  exactly  in  two  dimensions 
and  for  a  special  temperature  the  statistical  mechanics  of  the  interface 
between  two  charged  one  component  plasmas  separated  by  an  impermeable 
membrane. The  interesting  fact  about  this  model  is  that  it  tr'eats  the 
entire  double  layer , the*metal 0  and  the  *electroly te*  side, on 

an  equal  footing, and  it  can  be  switched  from  the  polarizable  to  nonpolar izable 
behaviour. The  above  calculations  could  be  easily  extended  to  more  complex 
situations. For  instance  in  Appendix  D,the  case  when  the  membrane  itself 
has  a  fixed  charge  density  is  discussed. 

Another  solvable  case  is  that  of  several  contiguous  impermeable 


membranes  separating  plasmas  of  different  charge  densities 


M 


Prom  the  microscopic  point  of  view, an  important  result  of  this  paper  is  that 
the  existence  of  an  impermeable  membrane  is  strictly  equivalent  to  a  step 
variation  in  the  fugacity  of  the  mobile  ions. For  the  permeable  membrane 
there  is  no  discontinuity,  and  therefore  the  potential  drop  across  the 
interface  is  constant .This  is  perhaps  the  first  rigorous  illustration 
of  the  fact  that  an  ideally  polarizable  electrode  exists  only  when  there  is 
no  miscibility  between  the  two  phases. 

In  conclusion  the  behaviour  of  this  model  exhibits  a  remarkable 
ressemblance  to  that  of  real  electrodes: the  potential  drop  across  the  interface 
is  the  natural  external  variable  (and  not  the  surface  charge  density^ 

Lippmann' s  equation  is  verified, the  shape  of  the  electrocapillarity 
and  capacity  curves  shares  common  features  with  experiment. 

The  remarkable  fact  that  such  widely  different  systems  have  similar 
characteristic  behaviour, seems  to  imply  that  they  are  not  very  sensitive  to 
the  structural  details  of  the  interface. 
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APPENDIX  A:  ELECT RONE JTRAL IT Y  CONDITION 


Consider  the  left  hand  side  of  the  interface  (x<0).From  (3.4) 


we  get 


/(A-a.A.  /  Jh  ^ 


Hence 
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Hence 
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Then,  using  (2.34)  we  readily  obtain  the  electroneutrality  relation  (3.10). 

The  procedure  is  analogous  for  x>0. 


APPENDIX  3s MULTI POLE  SUM  RULES 


The  first  and  second  multipole  sum  rules  are 


— oo  -4.0J  J 

J  d3  *1-3  V  ° 


Consider  first  (3.1) 
Prom  (3.15)  we  have 


where 


=  «.v/  W*!  /  dt 

llu,  U1  «oo  C/ 
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Hence 
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Performing  the  last  integration  we  find  (3.1). 
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Consider  now  (3.2)  ,which,  by  use  of  (B.l)  can  be  rewritten 
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We  have 
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The  last  term  between  brackets  is  the  derivative  of  the  denominator 
Integration  by  parts  finally  yields  (B.5). 
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APPENDIX  C: VERIFICATION  OF  SOM  ROLES 


We  first  consider  the  sun  rule 


from  (3.4)  we  get 


Jr  «V  [PC*}- £}]<** 

©  J 
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which  can  be  written 
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r=  -$—J±  /  dt- fc^c^mirVdA.ZA.V1}  -  AmhC  mh1  f  - C  --  — -  -  l]  O*) 

K*  V  *  Jmm9  <R  . 

Now  we  split  the  domain  of  integration  in  the  two  intervals 
[  -oo,  <r^  ],[  <r/v  ,  oo]  fand  we  regroup  the  terms  to  insure  convergence 
when  t->-oo.This  yields  ^ 

.si  r^' Jh r  -h’nht- 

1  L  *  J 

+  T  At  "icC"*)  ^  1  4  tV  (4  +  13^  (c-0 

Using  the  relations  (2*37)  ,(5. 3)  and  (5.4),  which  define  and 
the  saddle  point  equation  (2.34)  to  eliminate  , it  is  easy  to 

see  that  (C.4)  is  nothing  but  PT 

The  calculation  is  similar  for  the  left  hand  side  and  also  for  the 


symmetric  case  m-1 


Me  verify  now  contact  the  theorem  (5.18): 


Prom  (3.*)  we  have 
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which  can  be  written 
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After  similar  manipulations, ve  get 

OO 

=r  .  «Vj.  / _ c 

2Tt  Xv-  Jm0a  *  jr# 


(C-c) 
"J  <i  V- 


Bence  we  have 
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where  we  have  used 
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But  from  (3.4)  we  have  also 

TjVv  ^ 

Adding  (C.8)  and  (C.10)  we  finally  recover  the  contact  theorem  (5. It) 
fort4  -2  ^ 
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APPENDIX  D: CHARGED  MEMBRANE 

The  model  proposed  in  this  paper  can  be  extended  easily  to  the 
case  where  there  are  fixed  charge  distributions  on  each  side  of  the 
interface. 

Consider, for  instance, the  case  of  two  charged  rings, one  of  radius 

R  -  £  and  surface  charge  density  ,and  the  other  of  radius  Rw  +  £ 

*  'W 

and  surface  charge  density  €0^^, 

Let  £  ->0.This  is  equivalent  to  say  that  the  impermeable  membrane 

has  now  a  charge  density  «2L  »  -£LCG'7  +  0i  ^ 

7 1  7C 


The  calculation  follows  the  same  line  as  before  with  only  some 
minor  changes  (see  also  [11]) .The  electroneutrality  conditions  (2.1) 
and  (2.2)  become 


Xj-.  {*,*-«.* 


(3>.f  \ 


and  there  is  an  extra  term 

V  ’  *  -  «V  R  ZZ  'fn  C&-Z  \ 

“  i.w;.,  '•«.  1 

in  the  expression  (2.3)  of  the  potential  energy. 

In  (2.6)  we  now  have 

n?=  n  +  +  2(rR- 1 


and 


:  m.  ***  ‘  erfc  (><•-< 


instead  of  (2.13). 


Furthermore  we  now  have  icr  instead  of  G-  ,  in  all  the  integrals 

3.  * 

over  t.For  instance, the  saddle  point  equation  (2.34)  reads 


/ 


Ah  z 


is  now  a  function  of  Q*  and  0~  (and  not  only  of  0“  ) 


The  one  and  two  body  distribution  functions  are  still  given  respectively 
by  (3.4)  and  (3.15)  and  all  the  sum  rules  hold. 

Contrary  to  the  case  of  the  permeable  interface  (11), the  potential 
drop  across  the  interface,  depends  on  the  distribution  of  charges 

at  the  inter  face :the  surface  charge  distributions  of  the  interface  are  not 
screened  anymore. 
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FIGURE  CAPTIOUS 


Figure  lsDensity  profile  for  ra»0.5  and  Cfr-1. Units  are  e»*-*l.  (^s 
Figure  2:Level  curves  of  the  two  body  truncated  density  correlation 
function. Same  case  as  Figure  1. Particle  2  is  at  the  wall  on  the  RHS. 
Figure  3:3ame  as  figure  2. Particle  2  has  coordinates  (.5,0). 

Figure  4:Same  as  figure  2. Particle  2  has  coordinates  (1,0). 

Figure  5: Potential  drop  versus  surface  charge  density. The  primitive 
electrode  is  compared  to  the  full  coupled  model.  (m-.5) 

Figure  5 :Qif ferential  capacity  versus  potential  drop. Comparison  of  the 
primitive  electrode  to  the  full  coupled  model.  (m».  5) 

Figure  7 : Electrocap ilarity  curve,  (m*. 5)  .  (  p. *.c.  a  . 
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ABSTRACT 

A  microscopic  mode)  for  the  metal-solution  interface  is  presented,  which,  for  the  first  time,  includes 
both  a  non-ideal  treatment  of  the  metal  and  a  molecular  model  for  the  solution  side.  The  metal  is 
described  by  a  jellium  model  including  electron-ion  pseudo-potential.  The  solution  is  treated  as  a  mixture 
of  hard  spheres  with  point  dipoles  (solvent  molecules)  or  point  charges  (ions).  The  statistical  mechanics  is 
solved  using  mean  spherical  approximation.  No  a  priori  separation  is  introduced  between  compact  and 
diffuse  layers.  A  simplified  version  of  the  model  is  applied  to  the  case  of  Hg-DMSO  systems. 


INTRODUCTION 

The  classical  theories  of  the  ideal  polarized  electrode  are  based  on  two  premisses: 

(2)  The  metal  has  no  direct  effect  on  the  differential  capacitance  C,  although  it 
does  determine  the  potential  of  zero  charge  and  affects  the  orientation  of  solvent 
molecules.  From  an  electrostatic  viewpoint,  every  metal  acts  in  just  the  same  way, 
the  only  interaction  being  that  of  ideal  images. 

(2)  The  solution  side  has  two  layers:  the  external,  or  diffuse  layer,  which  can  be 
described  by  such  simple  models  as  the  Gouy- Chapman  (GC)  theory,  and  the 
internal,  or  compact  layer,  which  in  the  absence  of  specific  adsorption  is  formed 
only  by  solvent  molecules. 

These  two  premisses  are  now  being  challenged.  The  recent  important  progress  in 
the  theory  of  metals  and  metal  surfaces  makes  it  difficult  to  admit,  that  the  metal  is 
totally  unresponsive  to  changes  in  the  electrolyte.  The  surface  charge  distribution  in 
the  metal  must  respond  to  local  changes  in  the  electrostatic  Held  produced  by 
modifications  in  the  ionic  and  solvent  distributions.  This  will  surely  produce 
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contributions  to  the  differential  capacitance  that  will  be  different  for  each  metal, 
since  the  response  to  an  external  field  is  different  for  each  metal,  and  depends  on 
details  of  the  electronic  structure  near  the  surface.  In  previous  work  [1-3]  an 
estimate  of  these  effects  was  given.  The  solution  side  was  represented  by  a  layer  of 
dielectric  material.  In  a  recent  publication.  W.  Schmickler  [4]  has  discussed  a  similar 
model.  His  conclusions  about  the  role  of  the  metal  are  similar  to  ours. 

On  the  solution  side,  recent  advances  in  the  treatment  of  molecular  solutions  [5.6], 
i.e.  solutions  in  which  the  solvent  is  not  treated  as  a  continuum,  have  made  it 
possible  to  treat  the  statistical  mechanics  of  the  solution  side  of  the  electrode 
interface.  Although  the  approach,  in  principle,  can  handle  realistic  models  of 
solvents  and  ions,  the  initial  work  has  dealt  with  a  model  in  which  the  solvent  is 
represented  by  hard  spheres  with  point  dipoles  [7].  This  work  has  revealed  two 
important  deficiencies  of  the  classical  treatment: 

(1)  The  separation  of  the  inner  and  outer  layers  for  the  solvent  orientation  process 
(or  better  polarization)  implies  a  violation  of  Maxwell’s  equations.  The  dipoles  of  the 
solvent  are  oriented  wherever  there  is  an  electric  field,  i.e.,  throughout  the  entire 
electric  double  layer. 

(2)  The  use  of  local  dielectric  constants  is  also  unwarranted.  Although  this  is  a 
much  more  subtle  problem,  recent  theoretical  work  [8]  has  shown  that  dielectric 
constants  have  a  clear  meaning  only  for  large  systems  (thermodynamic  limit)  and 
not  at  the  molecular  level.  The  polarization  is  the  local  variable  that  is  required. 

In  this  work  we  present  a  comprehensive  theory  in  which  all  of  these  concepts  are 
included.  In  section  (I)  we  review  the  models  of  the  metal  and  solution*  sides  of  the 
interface.  The  potential  drop  and  capacitance  are  then  computed  using  a  simple 
GMSA  (generalized  mean  spherical  approximation)-based  model  which  has  the 
obvious  advantages  of  analyticity  and  simplicity. 

In  section  (II)  we  discuss  the  electronic  density  profile  of  this  model.  We  then 
consider  a  model  in  which  the  physical  parameters  are  those  of  an  alkali-ion  in 
DMSO  in  the  neighbourhood  of  a  mercury  electrode. 

In  our  model  the  point  of  zero  charge  is  not  the  point  of  zero  potential  (which  is 
the  case  in  the  restricted  primitive  model):  we  also  discuss  the  change  in  the  surface 
potential  of  the  metal  induced  by  the  presence  of  the  solution.  The  influence  of  some 
of  the  parameters,  and  detailed  comparison  between  Hg  and  Ga  electrodes,  will  be 
discussed  in  future  work. 

In  section  (V)  we  examine  possible  improvements  in  our  work. 

(I)  MODEL  FOR  SURFACE  INTERACTIONS 

Consider  a  plane,  ideally  smooth  electrode,  where  z  is  the  distance  from  the  metal 
surface.  The  bulk  metal  is  then  given  when  z  —  -  sc.  For  the  metal  we  use  a  jellium 
model  (with  inclusion  of  the  electron-ion  pseudo-potential).  The  details  are  given  in 
ref.  9,  but  for  completeness  we  will  give  an  outline  of  this  work.  The  metal  ions  are 
taken  as  a  continuum  of  density  p  for  z  <  0.  In  other  words 

p(z)-p0(-r) 
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where  9(x )  is  the  Heaviside  function 

jc)  —  1  x  >  0 

*(x)-0  x  <  0 

The  electronic  gas  is  free  to  move,  but  its  equilibrium  distribution  n(z)  will  be 
that  which  minimizes  the  total  surface  energy  of  the  system.  The  electron-ionic 
background  interaction  is  given  by  a  pseudo-potential.  The  properties  of  the  electron 
gas.  such  as  the  kinetic,  exchange  and  correlation  energies  are  computed  using  the 
local  approximation. 

On  the  solution  side  the  ions  are  represented  by  hard  spheres  of  diameter  a,  and 
charge  ±e  *,  and  bulk  density  p,.  The  solvent  is  also  represented  by  hard  spheres  of 
diameter  point  dipole  p  and  density  pd. 

We  assume  that  no  chemical  reaction  can  take  place,  and  no  miscibility  occurs 
between  the  phases  (10].  However,  this  model  does  not  exclude  the  existence  of  a 
surface  layer  in  which  ail  the  particles  are  present.  In  fact,  the  existence  of  this  layer 
is  consistent  with  contemporary  theories  of  the  metallic  surfaces.  Any  model  which 
does  not  include  this  layer  leads  to  predictions  which  are  incompatible  with 
experiments  (1.1 1J. 

In  our  model,  the  electrons  can  diffuse  into  the  ionic  solution.  In  actual  practice, 
calculations  show  that  both  in  the  interface  with  the  vacuum  [12],  or  with  a  dielectric 
film  [9],  the  electron  density  spill-over  is  very  small.  However,  the  change  in  the 
spill-over  with  the  electrode  charge  will  give  a  direct  contribution  of  the  metal  to  the 
differential  capacitance  [1,4]. 

The  metal  electrons  in  the  solution  are  subject  to  electrostatic  interactions 
(electron-ion  and  electron-dipole)  but  also  to  a  repulsion  due  to  the  core  electrons 
of  ions  and  solvent  molecules.  These  are  represented  by  the  Harrison  repulsive 
potential  (9]. 

There  are  several  parameters  in  our  model  that  are  fixed  a  priori,  but  which 
should  come  out  of  a  first-principles  calculation.  One  of  these  quantities  is  the 
distance  of  closest  approach  of  the  ions  and  solvent  to  the  plane  z  =  0  (the  metal 
surface).  We  will  assume  that  it  is  just  the  sum  of  the  molecular  radius  of  the  solvent 
or  ions  (which  for  this  simple  model  are  considered  to  be  equal)  and  the  radius  of 
the  metallic  ion.  In  reality,  this  parameter  should  arise  from  a  self-consistent  density 
functional  calculation  [13]. 

In  short,  we  will  assume  that  there  is  an  ideally  smooth  wall,  located  at  ;  =  5. 
which  cannot  be  penetrated  by  either  ions  or  solvent  molecules.  The  distances  of 
closest  approach  are 

‘,oa“5  +  «i/2  for  ions 

-  5  +  2  for  solvent  molecules 


*  We  use  Bohr  atomic  units:  t  -  electron  charge:  m  m  electron  mass;  e  »  m  -  1.  Bofh  radius:  a0  »  0.0519 
nm.  Energy  unit:  I  Hartree  «•  27.2  eV. 


:56 


There  is  also  the  problem  of  surface  smoothness  and  structure.  From  recent 
experimental  evidence,  it  is  known  that  certain  metal  surfaces  are  subject  to 
structural  changes  during  chemisorption  processes.  In  these  cases  the  surface  over- 
layer  changes  its  structure.  It  is  quite  conceivable  that  in  an  intrinsically  more  labile 
environment  such  as  that  of  liquid  metals,  the  metal  overlayer  would  also  respond 
structurally  to  changes  in  potential. 

Furthermore,  the  question  of  surface  ripples  (capillary  waves)  at  the  interface, 
and  its  effects  on  the  structure  and  thermodynamics  of  the  electrode  is  not  discussed. 
We  hope  to  address  these  and  other  questions  in  future  work. 

Summarizing,  the  ions  and  dipoles  in  the  solution  are  not  seen  simply  as  a  hard, 
charged  surface,  but  rather  the  sum  of  a  charged  surface  (located  at  5  +  a/2)  and  a 
distribution  of  electron  density  beyond  the  plane  of  closest  approach. 

Even  this  simple  model  is  not  easily  amenable  to  a  complete  numerical  calcula¬ 
tion.  We  will,  therefore,  in  this  first  communication,  use  a  simplified  theory  that  has 
the  advantage  of  giving  explicit  analytic  expressions  for  the  physically  relevant 
quantities. 

(II)  A  SIMPLIFIED  MODEL 

For  simplicity,  we  shall  discuss  only  the  restricted  model  in  which  n,  « <r*  « <r.  - 
<x4  ■  a  (although  not  in  published  form,  the  results  for  o  „  -  o  _  ~  <Jd  are  also  availa¬ 
ble).  We  also  neglect  the  effects  of  the  electron  density  on  molecules  and  ions  in  the 
solution.  We  shall  leave  the  discussion  of  the  validity  of  this  assumption  for  the 
future  (see  also  ref.  3  for  a  first  discussion  of  this  effect).  The  solution  is  represented 
by  the  MSA  [5,6]  which  is  only  qualitatively  correct.  Near  a  charged  interphase  a 
simplified  theory,  valid  presumably  for  dilute  solutions,  has  been  developed  by  Blum 


Fij.  1.  Simplified  model  of  the  interface:  4  -  position  of  the  ideally  smooth  wall;  a  -  diameter  of  the  ions 
and  solvent  molecules. 
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and  co-worker  (7J.  In  this  theory,  the  dipolar  interactions  are  decoupled  from  the 
ionic  distributions,  which  are  given  by  the  modified  Gouy-Chapman  (MGC)  theory. 
We  are  thus  satisfying  the  electrostatic  part  of  the  contact  density  theorem.  The 
simplified  model  of  the  interface  is  given  on  Fig.  1. 

The  charge  density  profile  is  given  by 

where  *  is  Debye's  inverse  length  and 
d,  m  S  +  c/2 

The  total  excess  charge  of  the  solution  side  is 

<?,-  f  dzqAz)-  -q m 

•v, 

The  electronic  density  profile  is  postulated  to  be  of  the  form  [1,9] 
n(r)-n[l 
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where  n  is  the  bulk  metallic  electron  density.  Clearly,  the  total  excess  charge  on  the 
metal  side  is 


<7m 


■ f  dz[p0(-z)-n(z)]- -nz0 

"  —  oo 


(4) 


This  last  relation  also  defines  r0. 

For  a  given  value  of  a,  and  ^M,  the  charge  and  polarization  (/>(r))  profiles  are 
computed  from  the  theory  [7],  These  in  turn  yield  the  potential  V(z)  and  electric 
field  £(z),  which  are  obtained  by  solving  Poisson’s  equation 


3£(z)  d2V{z)  ,  [a/  ,  ,  , 

){««,« . 


««0e  - - (5) 

The  potential  V(z)  is  continuous  in  all  space.  Because  of  the  discontinuities  in  the 
density  profiles,  V(z)  has  a  different  form  in  the  regions 

(I)  -ao  <  z  <  0 

(II)  0  <  z  <  d, 

(III)  d,  <  z  <  oo 

The  solution  of  eqn.  (5)  is  straightforward,  and  is  explained  in  Appendix  I.  The 
total  potential  drop  across  the  interface  is 


AK- K(oo)-  K-oo)-^ 
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The  integral  in  this  expression  is  explicitly  known.  The  reader  is  referred  to  ref.  7 
(eqn.  15)  for  the  derivation.  The  result  is 

•r  f  dzP(z)  -  lim  P(s) 

»-»o 


The  Laplace  transform  P(s)  is 
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Here  pd  is  the  bulk  number  density  of  solvent  molecules.  The  coefficients  y3h  and 
>3I2  are  functions  of  the  bulk  dielectric  constant  t  of  the  solvent,  and  are  given 
through  an  auxiliary  parameter  A 
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The  auxiliary  parameter  X  is  computed  from  the  equation 


t  »  AJ(1  +  A)4/16 


(8') 


The  functions  $,(r)  and  -4(r )  are  defined  by 

4>,(s)  -  l/(se)2[l  -  jo-e-”] 

^(j)  «  1  -24»i(5)(Xj-  1)-44.j(5)(Xj-X) 

with 
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Then 


Substituting  into  eqn.  (6)  yields 


(9) 
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at  the  point  of  zero  charge  (r0  -  0),  this  expression  simplifies  to 

-4 lrn/ci1  (12) 

The  potential  drop  across  the  interface  is  a  function  of  only  n  and  a,  which  are 
metal  parameters.  This,  however,  does  not  mean  that  the  potential  AK  is  equal  to  the 
surface  potential  of  the  metal-vacuum  interface  xm.  since  a  is  also  a  function  of  the 
solution  parameters.  In  the  limit  a  -*  oo  the  electronic  profile  tends  to  n(z)m 
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n9(z-zg),  and  then  eqn.  (11)  is  not  the  same  as  the  expression  given  in  ref  14. 
There  are  two  reasons  for  this  discrepancy.  First,  the  term  2itnzg  does  not  exist  in 
ref.  14,  because  now  the  charge  is  not  localized  at  the  plane  z  —  0.  Secondly,  the  third 
term  in  eqn.  (11)  is  not  that  of  ref.  14.  However,  if  d,  »■  a/2  then  they  are  identical. 
This  means  that  an  additional  capacitance  term  appears  because  of  the  shift  8  in  the 
metal  background  profile. 

The  differential  capacitance  of  the  system  is 

1  3AK  8AK  3A  V 

c"a*."  a?M“«ar0 


using  eqn.  (1 1),  we  obtain 

?“^[1+T(1+irl)]+4,r[5**01 
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In  this  expression  the  first  term  of  the  right-hand  side  is  the  result  of  the  MSA 
(7,8J.  If  the  solvent  molecules  are  shrunk  to  point  dipoles  ( od  -*  0),  one  should 
recover,  from  eqn.  (14),  the  differential  capacitance  of  the  MGC  theory.  Otherwise 
(14],  one  can  define  the  effective  dielectric  constant 


(15) 


which  then  yields  the  MGC  result  with  the  effective  dielectric  constant  i.  We  must 
remark  that  in  spite  of  the  formal  resemblance,  the  physics  is  different,  since  the 
solvent  is  polarized  throughout  the  entire  electric  double  layer. 

The  second  term  of  eqn.  (14)  arises  from  the  fact  that  the  electrode  charge  is  not 
localized  at  z  »  0.  The  next  term  is  due  to  changes  in  the  dipolar  surface  layer  of  the 
metal  as  a  function  of  the  excess  charge.  And  the  last  term  reflects  the  fact  that  the 
distance  d,  is  also  a  function  of  the  excess  charge.  A  detailed  discussion  of  this  fact 
can  be  found  in  ref.  3.  At  the  point  of  zero  charge  this  term  vanishes. 

Because  of  the  intrinsic  limitations  of  the  MSA,  we  will  limit  ourselves  to  small 
surface  charge  densities  around  the  point  of  zero  charge.  Thus,  for  a  given  value  of 
d(,  we  only  need  to  compute  a  for  the  differential  capacitance  C. 


(HI)  ELECTRONIC  DENSITY  PROFILE 


.  The  ansatz  (3)  on  the  electronic  density  profile  is  particularly  convenient  since 
n{z)  only  depends  on  the  parameter  «.  We  will  compute  this  parameter  by 
minimizing  the  surface  excess  energy  U,.  The  kinetic  energy  and  the  non-coulombic 
part  of  the  electron-neutralizing  background  are  not  affected  by  the  presence  of  the 
solution,  and  therefore  we  use  the  expressions  given  in  ref.  9  (eqn.  14  with  a  -  0). 
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Uk,n  +  U„e  +  Uinh  -  -[  -  1 .6423  *2/ 3  +  0.25037  «'/3  +  0.006584] 

+  a/i (log  2)/72  (16) 

The  pseudo-potential  contribution  Vpi  to  the  energy  is  given  by  (eqn.  10  of  ref.  1 
with  a  -  0) 


t/„, 


2ir/i: 
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where  <40  and  /?M  are  the  pseudo-potential  parameters  of  Heine  and  Animalu  (see 
also  ref.  1). 

The  electrostatic  energy  is  computed  from  the  charge  and  polarization  profiles, 
and  the  potential  V(z)  which  has  been  calculated  in  Appendix  I.  We  find 
.2 
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The  last  term  Q  consists  of  the  remaining  contributions  to  the  energy  which  are 
not  functions  of  a,  and  thus  irrelevant  to  our  calculation. 

The  penultimate  term  in  eqn.  (18)  represents  the  interaction  of  the  electrons  with 
the  local  polarization  of  the  solvent  molecules  in  the  solution.  This  contribution  is 
easily  calculated  from  the  Laplace  transform  P(s)  (eqn.  7): 

-111  _!l!e«o£(a)  (19) 

a  Jj  a 


There  is  also  another  contribution  due  to  the  repulsive  interactions  of  the 
electrons  and  the  electronic  clouds  of  the  ions  and  solvent  molecules  [9].  These 
interactions  are  represented  by  Harrison’s  pseudo-potentials,  with  parameters  X,  (for 
ions)  and  Xd  (for  solvent).  This  interaction  is  independent  of  the  ionic  charge  or 
solvent  dipolar  orientation,  and  because  of  the  use  of  the  MSA,  it  will  depend  only 
on  the  total  density  profile  pa(r).  Since  all  the  particles  are  hard  spheres  of  equal 
diameter,  a  single  parameter  X  (a  weighted  average  of  X,  and  Xd)  will  be  enough  for 
our  purposes.  We  obtain 


U%^\jirjiz'n{z)6{r-r’)^{z')  (20) 
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Ukin  +  +  £/inh--[  - 1.6423  n2/3  +  023037  n,/3  +  0.006584] 

CL 
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The  pseudo-potential  contribution  to  the  energy  is  given  by  (eqn.  10  of  ref.  I 
with  a  «■  0) 


+«(M  +  A,,i/1  -A,J>i/2  -  *„)  +  ^ 
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where  A0  and  are  the  pseudo-potential  parameters  of  Heine  and  Animalu  (see 
also  ref.  1). 

The  electrostatic  energy  is  computed  from  the  charge  and  polarization  profiles, 
and  the  potential  V(z)  which  has  been  calculated  in  Appendix  I.  We  find 
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The  last  term  Q  consists  of  the  remaining  contributions  to  the  energy  which  are 
not  functions  of  a,  and  thus  irrelevant  to  our  calculation. 

The  penultimate  term  in  eqn.  (18)  represents  the  interaction  of  the  electrons  with 
the  local  polarization  of  the  solvent  molecules  in  the  solution.  This  contribution  is 
easily  calculated  from  the  Laplace  transform  P(s)  (eqn.  7): 

r&tP(z)  #-*“*•>-  e"»^(a)  (19) 

cl  a 

There  is  »iy*  another  contribution  due  to  the  repulsive  interactions  of  the 
electrons  and  the'  electronic  clouds  of  the  ions  and  solvent  molecules  [9].  These 
interactions  are  represented  by  Harrison’s  pseudo-potentials,  with  parameters \  ,  (for 
ions)  and  Xtf  (for  solvent).  This  interaction  is  independent  of  the  ionic  charge  or 
solvent  dipolar  orientation,  and  because  of  the  use  of  the  MSA,  it  will  depend  only 
on  the  total  density  profile  p^(a).  Since  all  the  particles  are  hard  spheres  of  equal 
diameter,  a  single  parameter  X  (a  weighted  average  of  Xf  and  Xd)  will  be  enough  for 
our  purposes.  We  obtain 
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where  pd  is  the  dipole  concentration  and  p„  is  the  (  +  )  ion  concentration.  The 
integration  is  easily  performed: 

^8“Tjfd*e'a<'''0’po<2)  (21) 
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where  p0(a)  is  the  Laplace  transform  of  the  total  density  distribution  function 
(Po(*)“Pd(*)  +  P,(z)  +  p_(r)).  w*  obtain 
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The  functions  <p,(a)  and  <pj(a)  are  defined  by  eqn.  (9). 
The  total  surface  energy  is  then 

U,  -  «4,.  +  +  i/p.  +  i/.,  +  Up 

Using  eqns.  ( 1 6)— (22)  we  obtain,  finally. 
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The  coefficients  C,,  C2,  £  and  D,  are  given  in  Appendix  2. 

In  eqn.  (23),  all  the  electrostatic  interaction  terms  with  the  solution  cancel  at 
r0  •  0,  the  point  of  zero  charge.  The  shift  in  the  surface  potential  of  the  metal,  when 
it  is  taken  from  the  vacuum  to  the  solution  is  due  to  the  electron-solution  repulsive 
interactions. 

The  minimization  of  U,  then  yields  a.  This  parameter  is  then  used  to  compute  the 
properties  of  the  electrode  with  the  aid  of  eqns.  (11)  and  (19). 


TABLE  1 

Values  of  the  profile  parameter  a  (Hg)  as  a  function  of  <f,  and  X.  For  the  bare  surface  a0  -  0.901  a.u. 
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Two  of  the  parameters  of  the  problem  are  free,  and  have  to  be  selected  using 
physical  intuition.  They  are  S  and  X.  It  seems  realistic  to  choose  the  value  of  S 
between  0  and  Rc,  the  crystallographic  radius  of  the  ions  in  the  metal  {9). 

For  water  the  parameter^  has  been  estimated  to  be  -  15.  For  DMSO  we  might 
expect  an  even  higher  value  since  there  are  54  electrons  in  DMSO  against  10  in 
water.  In  Table  l  we  show  the  influence  of  X  on  a  assuming  that  -  <r/2  or 
dt  m  Re  +  o/l.  As  expected,  the  repulsion  term  with  the  solution  causes  an  increase 
in  the  values  of  a.  This  is  similar  to  what  we  found  previously  [9],  In  Fig.  3  we  have 

plotted  the  shift  in  the  surface  potential  of  the  metal  5xm**4irn(-^r — 7)  calcu- 

«  <*o 

lated  at  the  point  of  zero  charge  as  a  function  of  the  distance  d,.  The  corresponding 
case  of  gallium  for  \  -  30  is  also  shown. 

The  repulsion  from  the  solution  is  not  enough  to  keep  the  electrons  in  the  region 


Fig-  3.  Dependent*  of  txm  on  the  parameter  4,  for  different  choices  of  X:  (1)  Hg  X  - 15;  (2)  Hg  X  -  30; 
(3)  HgX'-4J;(4)OaX-30. 


3 


3 


di/au. 


Fig.  4.  Dependence  of  the  chugs  ipiU-over  if  on  the  parameter  d,  for  Hg  end  Ga  (X  -  30  and  c  -  0.1 
mot  I'1). 


z  <  dx.  There  is  always  a  small  charge  spill-over,  -  5 q,  in  the  region  z  >  dv 

-5<j-  -  /’**dr£e-« >*-*•»- 

'  Jd,  2  2a 

In  Fig.  4  we  show  iq(dt)  tor  both  Hg  and  Ga.  The  spill-over  8q  is  larger  for  Ga 
than  for  Hg. 

In  the  simple  model  presented  here,  this  effect  has  been  neglected  in  calculating 
the  charge  distribution  and  the  polarization  in  the  solution.  It  is  clear  that  this 
approximation  is  leu  justified  for  Ga  than  for  Hg.  However,  the  mathematical  tools 
to  include  it  are  available  and  we  will  indeed  discuss  these  effects  in  the  near  future. 

The  fact  that  there  are  electrons  on  the  solution  side  even  at  the  point  of  zero 
charge,  indicates  that  the  solvent  molecules  must  be  polarized  at  the  interface  even  at 
the  point  of  zero  charge.  Note  that  the  ratio  8q(Q«,)/tq(Hq)  is  not  simply  the  ratio 
of  the  electronic  densities. 

The  differential  capacitances  for  the  mercury  electrode  are  given  in  Table  2.  They 
were  computed  from  the  curve  bV m  f($M)  by  nume  Seal  differentiation.  With  the 


TABLE 2 

Differential  capacitance  of  the  interface  C/n F  cm'3  as  a  function  of  <f,  and  X.  The  concentration  it 
c  “  0.1  mol  1* 1 


X 

4, 

*«-! 

0 

j  +0.1  nm 

0 

2 

15 

9 

20.8 

100 

30 

9 

18.7 

58.8 

45 

8.7 

16.6 

43.4 

parameters  of  this  article  the  zero  charge  capacitance  corresponding  to  an  ideal 
metal  would  be 


4  ir  ,  ko 

l  +  T 


'  7.94  jiF  cm' 


We  see  that  the  total  capacitance  depends  crucially  on  the  distance  </, . 

When  d,  is  large  the  value  of  C  is  greater  than  C,  but  the  effect  of  non-ideality  of 
the  metal  remains  small.  On  the  other  hand,  when  the  distance  of  closest  approach 
decreases  (for  instance  if  there  is  interpenetration  between  the  metallic  ions  and  the 
particles  of  the  solution)  this  effect  becomes  quite  important  and  the  influence  of  X  is 
not  negligible  (for  comparison  the  experiment  value  of  18  /xF  cm-2  for  LiCIO,  ions 

I15D- 

First  calculations  show  that  the  differential  capacitance  for  Ga  would  be  greater 
than  for  Hg,  as  observed  experimentally. 


(V)  CONCLUSIONS 


We  have  presented  a  realistic  model  of  the  metal-solution  interface,  in  which  the 
metal  is  not  an  ideal  smooth  metal,  but  a  jellium  model.  A  simple  discussion  of  this 
model  is  given.  The  results  confirm  earlier  calculations  [1-4]  that  indicate  that  the 
metal  makes  an  important  contribution  to  the  differential  capacitance  C. 

The  solution  side  of  our  model  is  a  mixture  of  hard  spheres  with  charges  and 
point  dipoles.  No  a  priori  separation  of  compact  and  diffuse  layers  is  postulated  in 
this  model,  and  it  also  does  not  involve  the  use  of  local  dielectric  constants  (a 
meaningless  quantity).  The  calculations  show  that  the  solvent  molecules  are  polarized 
even  at  the  point  of  zero  charge.  This  is  so  because  of  the  spill-over  of  electronic 
charge  into  the  solution,  which  produces  non-zero  electric  fields  in  the  metal-solu¬ 
tion  boundary  layer.  Although  we  have  treated  this  effect  only  crudely,  the  recent 
work  on  the  GMSA  [16]  of  the  ion-dipole  mixture  allows  a  consistent  calculation  of 
this  effect 

We  also  show  that  the  distance  <f,  of  closest  approach  of  the  metal  ions  and  the 
solution  ions  is  a  crucial  parameter  in  the  determination  of  the  differential  capaci¬ 
tance  C.  However,  this  is  also  a  drawback  since  the  fact  that  d,  is  unknown  makes 


the  comparison  of  C  with  experiment  a  not  very  meaningful  exercise.  Indeed  a  value 
of  d,  which  fits  the  experiments  could  be  found. 

However,  the  model  that  we  have  discussed  here  represents  the  most  comprehen¬ 
sive  one  that  we  are  aware  of.  We  are  undertaking  a  detailed  study  of  the  influence 
of  various  parameters,  such  as  ionic  concentration,  density  and  dipole  moment  of 
the  solvent,  etc.  on  the  differential  capacitance.  The  results  will  be  published  in  the 
near  future. 
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APPENDIX  1 

Electrostatic  potential  profile 

Integrating  eqn.  (S),  obtain  (also  using  eqns.  1  and  3) 

(a)  For  negative  metal  charge  qu  (z0  >  0): 

V{z)—  -  (-eo<z<0) 

a 1 

V(z)  -  2irn(zi  -  t'+’—'/a2)  (0  <  z  <  z0) 

F(z) «  2irn[z0(2z  —  z0)  —  2/a2  +  e~*<‘~,,>/a2]  (z0<z  <d,) 

V(z)-2vn\z(t{2di  -r0)-2/a2  +  e_,K,_**)/«2]  +4 

+4*f‘dz’P(z’)-tm±e-<‘-‘0  ( z>dt ) 

Je,  * 

(b)  For  positive  qM  (z0  <  0): 

F(r)- -•—e^*— »  (-»<z<r0) 

a 

V(z)  m  -2»n[(z -z0)2  +  2/a2  -  e~*<,“,,,/a2]  (z0  <  z  <  0) 

F(z)  -  2wn[z0(2z  -  z0)  -  2/a2  +  e"«<‘—>/a2]  (0  <  z  <  d,) 

For  d,  <  z  <  oo  the  equations  are  those  of  the  case  q M  <  0. 

APPENDIX  : 

We  give  explicit  formulae  for  the  coefficients  of  eqn.  (23): 

C,  -  —  1 .6423n2/I  +  0.25037* 1/2  +  0.006584] 
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1  log  2 
irn  72 
X  (l+2it) 

W2Pd 

m(2  +  X) 

These  coefficients  are  computed  once  and  for  all  for  a  given  metal  and  solvent. 
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